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Abstract 

We study, by numerical simulations on a lattice, the behaviour of the gauge-invariant 
quark-antiquark nonlocal condensates in the QCD vacuum with dynamical fermions. A 
determination is also done in the quenched approximation and the results are compared 
with the full-QCD case. The fermionic correlation length is extracted and compared with 
the analogous gluonic quantity. 
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1. Introduction 



The so-called "nonlocal condensates", i.e., gauge-invariant field correlators, are the orig- 
inal starting points for any standard calculation which adopts the method of QCD sum 
rules Q . These quantities appear when evaluating the power corrections, via the "Oper- 
ator Product Expansion" (OPE) 0, of the product of two hadronic currents. The effects 
due to the x-distribution of these vacuum fluctuations sometimes have been neglected, 
dealing only with local condensates (g(O)g(O)), (G(O)G(O)), etc.. However, it has been rec- 
ognized in 0, £|, ||, H [7|, H U that in many applications the effects due to the ^-distribution 
of the nonlocal condensates have physical relevance and cannot be neglected. Therefore, 
the knowledge of those nonlocal condensates from first principles can be important for 
the study of the strong interaction theory and its applications. 

In a previous series of works, we have studied, by numerical simulations on a lattice, 
the gauge-invariant two-point correlators of the gauge field strengths in the QCD vacuum: 

2W(z) = (0\Tr{G lip (0)S(Q,x)G u<T (x)^(Q,x)} |0> . (1.1) 



G> p = 9T a G a w is the gauge field-strength tensor and 5(0, x) is the Schwinger phase 
operator needed to parallel-transport the tensor G va (x) to the point 0. 

These correlators have been determined on the lattice in the quenched (i.e., pure 
gauge) theory, with gauge-group 577(2) [fLOR , in the quenched 577(3) theory in the range 
of physical distances between 0.1 and 1 fm [11| [12| and also in full QCD, i.e., including 



the effects of dynamical fermions Jug. The basic results of all these determinations is 
that the correlator "D^p^^x), in the Euclidean theory, can be written as the sum of a 
perturbative-like term, behaving as l/|x| 4 , and a nonperturbative part, which falls down 
exponentially 

2^(*)~exp(-|*|/A A ) . (1.2) 

The correlation length is A^ ~ 0.13 fm for the SU{2) pure-gauge theory |T0| , A^ ~ 0.22 
fm for the SU(3) pure-gauge theory [jll, |12| and A^ ~ 0.34 fm for full QCD (approaching 



the chiral limit) |g. 

Along the same line, in this paper we present a lattice determination of the quark- 
antiquark nonlocal condensates (g(0)5(0, x)q(x)), adopting the same basic strategies and 
techniques already developed for the study of the gluonic correlators (|TJ, ITJL O. In 
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particular, we shall make use of the "cooling" technique [14], |15| in order to remove 
the effects of short-range fluctuations on large-distance correlators and get rid of the 
renormalizations. We shall not present here again the details of our cooling procedure, 
for which we refer the reader to our previous works O, [T2|, f 
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The quark-antiquark nonlocal condensates have been determined both in the SU(3) 
pure gauge theory and in full QCD, i.e., including the effects of dynamical fermions. The 
details of the computations and the results are presented in Section 2. In Section 3 we 
conclude with some remarks about the results. 



2. Computations and results 

As we shall explain in more detail later, for our lattice computations we have used four 
flavours of staggered fermions, so that we are considering a theory with four degenerate 
quark flavours in the continuum limit. Therefore, we have decided to consider only the 
expectation values of those quark-antiquark operators which are diagonal in flavour, but 
are nontrivial with respect of the Dirac spinor indices, i.e.: 

Q(x) = -j2(Tr[qm(r) ab S(0,x)q f b (x)]) . (2.1) 
/=i 

A few words about the notation used in (2.1). S(x, y) is the Schwinger string from x to y 

S(x,y) = Pexp(i9 j" dz» A^z)} , (2.2) 

needed to make Ci(x) a gauge-invariant object. "P" stands for "path ordering": for 
simplicity, we take S(x, y) along the straight-line path from x to y. = A a ^T a is the 
gauge field operator and T a are the matrices of the algebra of the colour group SU(N C ) 
in the fundamental representation (For N c = 3, T a = X a /2, where A a are the Gell-Mann 
matrices). The trace in (2.1) is taken with respect of the colour indices. The index / 
in (2.1) is a flavour index (/ = 1,2,3,4). The matrices P are the sixteen independent 
4x4 matrices of the Clifford's algebra acting on the Dirac indices a, b: 1, 7^, 7^, 7^7^, 
[7e,7e]/2, where 7^ (/x = 1,2,3,4) and 7^ are the Euclidean Dirac matrices. 
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Making use of the invariance of the theory under P, T, translations and rotations, 
one can easily verify that most of the sixteen vacuum expectations values (2.1) are zero 
for symmetry reasons. One is left with only two nontrivial quark-antiquark nonlocal 
condensates, out of the sixteen quantities Cj(x) in (2.1), namely 

Co(\x\) =-J2(Tr[Q f a(Q)S(0,x)q f a (x)]) ; 



C v (\x\) = ^C M (x) = -gy E(^(0)(«M)9!(i)]) • (2-3) 
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For simplicity, in the following we shall refer to these two quantities as, respectively, the 
"scalar nonlocal condensate" [C (\x\)] and the "(longitudinal-) vector nonlocal conden- 
sate" [C v (\x\)\. 

In order to construct our operators on the lattice we have used the following proce- 
dure. In the staggered formulation quark fields live on elementary hypercubes, so our 
correlators can be defined only for an even distance d in lattice spacing units. In com- 
puting (g(O)S'(O, x)q(x)) we have always put the point in the hypercube at the origin of 
the lattice, while x has been varied along the coordinate axes. The staggered propagators 
(XiXj) have been computed at first and have been connected to the origins of the hyper- 
cubes by minimal paths of gauge links (an average over paths of equal length has been 
performed). They have then been combined to build up the quark propagator. Finally, 
the Schwinger line connecting the origins of the two hypercubes has been put in. We have 
also performed an average over different directions. 

In this way we can construct two (adimensional) lattice operators C^d-a) and C^(d-a) 
(x = d ■ a, where a is the lattice spacing and d is the number of lattice spacings), which 
are proportional to Co(d ■ a) and C v (d ■ a) respectively in the naive continuum limit, i.e., 
when the lattice spacing a — > 0: 

C^{d ■ a) ~ a 3 C {d ■ a) + C(a 4 ) , 

CHd ■ a) ~ a 3 CJd ■ a) + 0(a 4 ) . (2.4) 

a— >0 

Higher orders in a in (2.4) as well as possible multiplicative renormalizations are removed 
by cooling the quantum fluctuations at the scale of the lattice spacing, as explained in 



Refs. fLl , [T|, 13J, [14|, This removal will show up as a plateau in the dependence of 
the correlators on the number of steps of the cooling procedure: our data are the values 
of the correlators at the plateaux. 
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The computations have been performed with four flavours of staggered fermions and 
the SU(3) Wilson action for the pure-gauge sector: we have considered both the full 
QCD case (i.e., including the effects of dynamical fermions) and the quenched case (where 
the effects coming from loops of dynamical fermions are neglected, i.e., det-fTLA] = 1 in 
this approximation, K[A] being the fermions matrix). 

For the case with dynamical fermions we have measured the nonlocal condensates on a 
16 3 x 24 lattice at f3 — 5.35 (j3 = Q/g 2 , where g is the coupling constant) and two different 
values of the dynamical quark mass: a ■ m q = 0.01 and a ■ m q = 0.02. In both cases the 
quark mass used in computing the quark propagator has been chosen to be the same one 
used during the simulation. Our samples were ~ 300 configurations at a ■ m q = 0.01, 
each one separated by 9 molecular dynamics time units, and ~ 70 at a ■ m q = 0.02, each 
one separated by 12 molecular dynamics time units. The computation was done on an 
QH4-APE machine by a standard hybrid Monte Carlo algorithm. 

For the quenched case the measurement has been performed on a 16 4 lattice at /3 — 
6.00, using a quark mass a-m q = 0.01 for constructing the external-field quark propagator, 
and also at /3 — 5.91, using a quark mass a ■ m q = 0.02. In both cases the value of f3 was 
chosen in order to have the same physical scale as in full QCD at the corresponding quark 
masses, thus allowing a direct comparison between the quenched and the full theory. In 
this way we can see if the inclusion of dynamical fermions has or has not considerable 
effects on the quantities that we are measuring. 

In the quenchedca.se two other measurements have been done at f3 — 6.00, using quark 
masses a • m q = 0.05 and 0.10 for constructing the external-field quark propagator, in 
order to study the dependence of the nonlocal condensates on the valence quark mass. 

Finally, in all cases, we have measured the nonlocal condensates at distances d = 
2, 4, 6, 8 in units of lattice spacing. 

In Figs. 1 and 2 we display the results for a 3 Co(d ■ a) and a 3 C v (d ■ a) respectively, 
versus the distance d in lattice spacings, for full QCD at a ■ m q = 0.01 and for full QCD 
at a ■ m q = 0.02. The results for a 3 Co and a 3 C v obtained in the quenched determinations 
are displayed in Figs. 3 and 4, for the cases f3 = 6.00 and a ■ m q = 0.01, (3 = 5.91 and 
a ■ m q = 0.02, [3 = 6.00 and a ■ m q = 0.05, (5 = 6.00 and a ■ m q = 0.10. 

For the scalar nonlocal condensate we have tried a best fit to the data with the function 




(2.5) 
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Results obtained in the various cases are shown in Table I. 

Analogously, for the longitudinal-vector nonlocal condensate we have tried a best fit 
to the data with the function 



C v (x) = A v x 3 exp(— ji v x) + 




(2.6) 



Results obtained in the various cases are shown in Table II. 

(The values of x 2 /Nd.o.f. reported in Tables I and II should be considered as purely in- 
dicative of the goodness of our best fits, since Nd.o.f. = ndata—n pa ram. is only an upper limit 
to the effective number of degrees of freedom after taking into account the correlations 
between data at different distances.) 

The form of the perturbative-like terms in Eqs. (2.5) and (2.6) (i.e., B /x 2 for the 
scalar condensate and B v /x 3 for the vector condensate) is that obtained in the leading 
order in perturbation theory, in the chiral limit m q — > 0: 



for a theory with Nf flavours and N c colours. In our case Nf = 4 and N c = 3: the values 
for the coefficients aB and B v obtained in the best fits agree as an order of magnitude 
with the estimates (2.7). 

To conclude this section, let us make some remarks about the significance of the results 
obtained for the longitudinal vector nonlocal condensate. As the number of cooling steps 
tends to oo , the gauge-field configuration is driven towards the zero-field configuration, 
so we expect both Cq{<1) and C„(d) to reach asymptotically their free-field values. In 
both cases, and unlike the case of the gluon field-strength correlators, these free-field 
values are different from zero and correspond to the leading order in lattice perturbation 
theory. 

We have determined the free-field values by measuring Cq{<£) and C^(d) on the zero- 
field configuration. Results are displayed in Figs. 5 and 6 respectively and compared to 
results obtained in full QCD for a ■ m q = 0.01. Qualitatively similar figures are obtained 
for other values of a ■ m q and (5. 




(2.7) 
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It clearly appears that the vector nonlocal condensate is strongly dominated by the 
free-field signal: our method to determine it is thus likely to be poorly sensitive to the 
nonperturbative signal. Moreover, the free-field signal for the vector nonlocal condensate 
is only weakly mass-dependent (it is different from zero in the chiral limit): this fact, to- 
gether with the above-mentioned dominance of the free-field signal, gives an explanation 
of why the data of the vector nonlocal condensate, reported in Figs. 2 and 4, depend so 
weakly on the value of the (valence) quark mass. We regard results concerning the vector 
nonlocal condensate as preliminary: it will be the subject of further investigations. 

This is not the case of the scalar correlator: there the free-field signal is only a 
small fraction of the whole signal at most distances. Moreover, as a further check of the 
reliability of the results obtained in this case, we have tried to subtract the free-field 
signal from the measured scalar nonlocal condensate: we have obtained that after the 
subtraction the results of the fits are unchanged for the nonperturbative part, while they 
give a value compatible with zero for the perturbative-like term. 



3. Discussion 



In principle from our simulations we can extract two quantities of physical interest: the 
quark condensate (q(0)q(0)) and the correlation length Ao = l///o of the scalar quark- 
antiquark nonlocal condensate. As explained in Refs. ||, j[| §, [7|, |], [| , Ao plays a relevant 
role in many applications of QCD sum rules, especially for studying the pion form factors 
and the pion wave functions. 

From the lattice we have extracted Ao in units of the lattice spacing a. To convert 
from these units to physical units, the scale must be set by comparison with some physical 
quantity. This is usually done by computing the string tension and/or the p mass on the 
lattice and comparing them with the physical values. In the full-QCD case, on the 16 3 x24 
lattice at /3 — 5.35 with four flavours of staggered fermions, we have found, by measuring 
m n and m p on our configurations and following the same procedure described in Ref . [ 16| , 
the following values for the lattice spacing: 

a(/5 = 5.35) ~ 0.101 fm, for a ■ m q = 0.01 ; 

a(/3 = 5.35) ~ 0.120 fm, for a ■ m q = 0.02 . (3.1) 
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In the quenched case the lattice spacing is approximately \\[% : 

a {YM \(3 = 6.00) ~ 0.103 f m ; 

a (YM \(3 = 5.91) ~ 0.120 fm . (3.2) 

Using these values, we can extract the physical values of the correlation length Ao for all 
the cases examined, obtaining the results reported in Table III. 

At a-rriq = 0.01 the value is roughly twice as big as the value for the correlation length 
A^ of the gluon field strength, defined in Eq. (1.2), at the same quark mass a ■ m q = 0.01 
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The full-QCD and the quenched correlation lengths Ao are nearly the same, when 
compared at the same quark mass. 



Similarly to the gluon case |Tj|, the fermionic correlation length appears to decrease 
when increasing the quark mass. However, our results seem to suggest that the value 
of A is sensitive to the value of the valence quark mass used in constructing the quark 
propagator, but is not much influenced by the dynamical quark loops coming from the 
determinant of the fermionic matrix det(-fT[A]). 

Using the values of the pion mass m n , measured on our configurations in the full-QCD 
case [a ■ m n = 0.263(8) for a ■ m q = 0.01, and a ■ m n = 0.372(13) for a ■ m q = 0.02], and 
the corresponding values of a ■ //o reported in Table I, we find that m n /fi = 1.6(4) for 
a ■ m q = 0.01 and m n /fio = 1.4(2) for a ■ m q = 0.02. Therefore, the inverse of the scalar 
correlation length, fiQ = 1/Ao, turns out to be proportional (within the errors) to the pion 
mass m^. in other words, it is the pion mass m v which determines the length-scale of the 
scalar nonlocal condensate. 

We conclude with a brief comment about the quark condensate. A way to extract this 
quantity is to consider the uncooled values of the scalar quark correlator at zero distance, 
Cq(0), for different quark masses and then to extrapolate those values to zero quark mass. 
From our full-QCD simulations we have obtained the following results (in lattice units): 

C L (0) = 7.17(6) x 10~ 2 , for a ■ m q = 0.01 ; 

C L (0) = 1.26(2) x 10"\ for a ■ m q = 0.02 . (3.3) 

A linear extrapolation of these two values to zero quark mass leads to the result: 

C L (0)U 9=0 = 0.0174(32) . (3.4) 
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Adopting the same procedure outlined in Ref. |16| (which properly takes into account the 
anomalous dimension of the quark condensate), we can then extract the following value 
for the single-flavour quark condensate in the MS renormalization scheme at a scale of 
H = 1 GeV: 

\(uu)\ m (fi = 1 GeV) = 0.013(2) GeV 3 . (3.5) 
This results is in perfect agreement with the phenomenological value |18 



As an alternative method, one could also try to extract the quark condensate directly 
from the values of the coefficient A , defined in Eq. (2.5). However, the values of a 3 A 
reported in Table I in the full-QCD case for a-m q = 0.01 and a-m q = 0.02 are visibly much 
smaller than the corresponding values for Cq (0) reported in Eq. (3.3). Apparently, no 
reasonable quark condensate can be extracted from these values of a 3 A . The reason for 
this discrepancy could lie in the anomalous dimension of the quark condensate: while we 
believe that the correlation length Ao is not affected by the cooling procedure, we do not 
know the effects of cooling on the quark condensate, which has an anomalous dimension. 
(This problem is not present in the case of gluon correlators, since the extracted gluon 
condensate is renormalization group invariant.) We hope to come back to this point in 
future works. 
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TABLE CAPTIONS 



Tab. I. Results obtained from a best fit to the data of the scalar nonlocal condensate with 
the function (2.5), in the various cases that we have examined ("/" stands for "full- 
QCD", while "g" stands for "quenched"). 

Tab. II. Results obtained from a best fit to the data of the longitudinal-vector nonlocal 
condensate with the function (2.6), in the various cases that we have examined ("/" 
stands for "full-QCD", while "g" stands for "quenched") . 

Tab. III. The physical values of the correlation length Ao for all the cases that we have ex- 
amined. Reported errors refer only to our determination and do not include the 
uncertainty on the physical scale. 
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Table I 



f3, theory 


a ■ m q 


a 3 A x 10 2 


a/j, 


aB x 10 1 


X 2 /N d . .f. 


5.35, f 


0.01 


0.49(13) 


0.16(4) 


0.13(3) 


1.3 • 10- 2 


5.35, f 


0.02 


1.7(5) 


0.26(4) 


0.19(10) 


6.4- 10" 3 


6.00, q 


0.01 


1.6(5) 


0.16(4) 


0.09(12) 


7.6 • 10" 2 


5.91, q 


0.02 


2.3(7) 


0.26(3) 


0.25(14) 


5.2 • 10" 2 


6.00, q 


0.05 


1.8(4) 


0.34(2) 


0.7(1) 


0.2 


6.00, q 


0.10 


5.6(5) 


0.55(1) 


1.0(1) 


1.3 - 10" 2 



Table II 



(3, theory 


a ■ m q 


a 6 A, 


a/i v 


B v 


X 2 /N d . . f . 


5.35, f 


0.01 


0.19(1) 


1.46(1) 


2.87(3) 


1.76 


5.35, f 


0.02 


0.20(1) 


1.48(2) 


2.84(3) 


0.8 


6.00, q 


0.01 


0.21(1) 


1.48(1) 


2.83(2) 


0.4 


5.91, q 


0.02 


0.216(7) 


1.494(6) 


2.81(1) 


2.3- 10" 3 


6.00, q 


0.05 


0.22(1) 


1.50(1) 


2.80(1) 


2 • 10~ 2 


6.00, q 


0.10 


0.26(1) 


1.54(1) 


2.73(2) 


2.2 



Table III 



P, theory 


a ■ m q 


A (fm) 


5.35, f 


0.01 


0.63^ 


5.35, f 


0.02 


0.46^ 


6.00, q 


0.01 


0.64±U; 22 


5.91, q 


0.02 


u - 4O -0.05 


6.00, q 


0.05 


0.30(2) 


6.00, q 


0.10 


0.187(3) 
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FIGURE CAPTIONS 



Fig. 1. The function a 3 Co(x) (scalar correlator) versus the distance d = x/a in lattice 
spacings, for the full-QCD case at f3 — 5.35 and quark masses a-m q = 0.01 (circles) 
and a • m q = 0.02 (squares). The curves correspond to our best fits [Eq. (2.5)]. 

Fig. 2. The function a 3 C v (x) (vector correlator) versus the distance d — x/a in lattice 
spacings, for the full-QCD case at f3 — 5.35 and quark masses a-m q = 0.01 (circles) 
and a ■ m q = 0.02 (squares). This last set of symbols (squares) has been shifted 
horizontally in the right direction to distinguish it from the other set of symbols 
(circles). As an example, we plot the curve corresponding to our best fit [Eq. (2.6)] 
to the data at a ■ m q = 0.01. 

Fig. 3. The same as in Fig. 1 for the quenched case at (3 — 6.00 and a quark mass a - m q = 
0.01 (circles), (3 = 5.91 and a ■ m q = 0.02 (squares), (3 = 6.00 and a ■ m q = 0.05 
(triangles down), (3 = 6.00 and a ■ m q = 0.10 (triangles up). 

Fig. 4. The same as in Fig. 2 for the quenched case at (3 — 6.00 and a quark mass a - m q = 
0.01 (circles), (3 = 5.91 and a ■ m q = 0.02 (squares), (3 = 6.00 and a ■ m q = 0.05 
(triangles down), f3 = 6.00 and a ■ m q = 0.10 (triangles up). Again, different sets 
of symbols have been shifted horizontally in the right direction. As an example, 
we plot the curve corresponding to our best fit [Eq. (2.6)] to the first set of data 
(circles). 

Fig. 5. The free-field value of Cq (d), measured on the zero-field configuration (filled cir- 
cles), compared to the results obtained in full QCD at a ■ m q = 0.01 (open circles). 

Fig. 6. The free-field value of C^(d), measured on the zero-field configuration (filled cir- 
cles), compared to the results obtained in full QCD at a ■ m q = 0.01 (open circles). 
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